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Abstract 

We investigate the exact analytic solutions for the growths of the 
dark matter and the baryon in sub-horizon scale. The growth of the 
dark matter 5dm is related to that of the halos. Thus, the exact so- 
lution for the growth of the dark matter is important to obtain the 
proper properties of dark matter halos. However, the dark energy 
model dependence of 6dm is confused with the (5dm dependence on 
fij^. Thus, the careful investigation is necessary for the 5dm depen- 
dence on dark energy models. We also obtain the exact solution of the 
growth of the baryon 5b which can be used to obtain the baryon bias 
factor h(a). This might be able to be observed in intracluster gas or 
in Lyman-a clouds. However, b{a) is quite model independent. Re- 
cently, we obtained the exact analytic solution for the growing mode 
solution of the matter linear density perturbation 5 on sub-horizon 
scale for general dark energy model [5]. This solution is not same as 
the well known approximate analytic solution [7]. The exact analytic 
solution shows the same evolution behavior of the growth factor ob- 
tained numerically. However, the exact solution is simple and useful 
for the extension to other models including modified gravity theories. 
Furthermore, it guides to the fact that the growth index parameter 
depends on both models lo^c and and thus we need to be careful 
for applying the fitting formulae to the general models [10]. The exact 
analytic solutions for the growth factor will provide the more accurate 
tools for the weak lensing, the number density of clusters, their mass 
and etc. 
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In the concordance model of cosmology baryon contributes about 4% of 
the total energy density with 23% cold dark matter and 73% dark energy 
in the universe [1]. The linear perturbation of the dark matter is related 
to the halos density perturbation. Thus, the exact solution for the dark 
matter growth Sdm provides the accurate and powerful tool for the probing 
dark matter halo properties. Also after the decoupling from the photon, the 
linear perturbation of the baryon Sb grows in the gravitational potentials 
created by the dark matter. Zero pressure of dark matter causes the uniform 
growth of the dark matter and the baryon in sub-horizon scale. Thus, the 
growth of baryon is related to the growth of the dark matter through the 
baryon bias factor Sb = b{a)dr)M [2]. The baryon bias is able to be observed 
when we consider not only the baryon distribution clustered in luminous 
objects but also the baryon not collapsed in galaxies. 

We use the flat Priedmann-Robertson- Walker universe to probe the linear 
density perturbations (5dm and 6b 
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where LOde is the equation of state (eos) of dark energy, pcr is the critical 
energy density, pnu, Pb and pde are the energy densities of the dark matter, 
the baryon and the dark energy, respectively. We consider the constant Ude- 
The sub-horizon scales linear perturbation equations for the coupled system 
of the dark matter and the baryon with respect to the scale factor a is given 
by [2, 3] 
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where we use the fact that the dominant contribution to the gravitation 
potential comes from the dark matter (i.e. J^dm<5dm + ^b^b — ^^dm^m)- 
We only consider the matter perturbation in Poisson equation because the 
dark energy is supposed to be smooth on sub- horizon scale [4] . Recently, we 
obtain the exact analytic solution of sub-horizon scale growth factor of the 
matter with the general dark energy models [5]. Compared to the previous 
results in Ref. [5], we just consider the dark matter as the source of the 
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gravitational potential. Thus, the solution for the sub-horizon scale linear 
perturbation of the dark matter Sdm is almost identical with that of total 
matter given in the reference [5] 
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where F = ^a^^^- = Qa^"^^-, = O?,^ + Q% and R = We use 

i? — > 1 (i.e. rij^jyj — >■ Q'^) in the second approximation. Thus, the second 
equality is valid as long as the dark matter energy density dominates the 
baryon energy density. We use ^ij^j^ = 0.23, fig = 0.4 and $7°^ = 0.73 later 
to be consistent with the concordance model. The baryon contribution is 
still about 17% of the dark matter. Thus, we use the solution of Jdm given 
in Eq (5). As in our previous work [5], we need to determine the coefficients 
of (5dm to make the above analytic solution as a growing mode one. We 
obtain them from initial conditions of growing mode solution 

<5DM(ai) ~ ttj and =i. (7) 

da ai 



We show the behaviors of the DM growing mode solutions in Fig. 1. 
In Fig. la, we show the evolutions of growing mode solutions of 5dm for 
the different dark energy models when 0dm = 0.23. The dotted, solid and 
dashed lines correspond to Ude = —1.5, —1.0 and —0.8, respectively. If there 
is the same amount of the dark matter at the present epoch for the different 
dark models, then there will be more matter component in the past for the 
smaller value of Ude- Thus, dark energy model with the smaller value of 
ojde maintains the longer linear growth behavior of (5dm- We show the O^j^ 
dependence of (5dm for the ACDM model (i.e. with cOde = —1.0) in Fig. lb. 
The dotted, solid and dashed lines correspond to 0,^)^,^ = 0.26, 0.23 and 0.20, 
respectively. Again if there is more dark matter at present epoch for the 
same models, then it maintains the longer linear growth behavior. The 5dm 
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Figure 1: a) Evolutions of 5dm(o) when = 0.23 for uj^e = —1.5,-1.0, 
and —0.8, respectively (from top to bottom), b) (5dm (a) evolutions for the 
different values of = 0.26, 0.23, and 0.20, respectively (from top to 

bottom). 



behaviors for the different dark models in Figs la look quite similar to (5dm 
of the cosmological constant model with the different values of in lb. 
It means that we have ambiguities in probing dark energy models with only 
using the growth of the DM. 

In addition to this, we are able to find the exact analytic solution for the 
growth of the linear perturbation of baryon 6-q from the above equation (4). 
The solution of 5-q consists of two components because the evolution equation 
of (5b is an inhomogeneous second order differential equation. Thus, we can 
separate the solution into the homogeneous solution and the inhomogeneous 
(particular) one 

(5b = (5^ + , (8) 
da'^ \ da a I da ' 

The equation of the inhomogeneous solution (5g is identical with that of (5dm 
in Eq. (3). Thus, the inhomogeneous solution of 6b is nothing but (5dm. We 
only need to solve the homogenous equation to obtain a complete solution 
of (5b. (5b is obtained from the similar method used in the previous works 
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where Chi and Ch2 are the integral constants. However, we are not inter- 
ested in the constant solution and remove Ch2 from the solutions. Thus the 
complete set of linear perturbation of the baryon is given by 
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where h{a) is the baryon bias factor. Again, the analytic solution of (5b (a) 
in Eq. (12) does not have any physical meaning before we fix the coefficient 
Chi- Thus, we use the fact that (5b — at decoupling epoch Odec ~ 0.001. 
Actually, we use Odec = 0.01 in order to include more relevant modes and 
more accurate to the matter dominated epoch approximation. Then we are 
able to obtain Chi from 

'^B(adec) = <^B(«dec) + '^DM(adec) = . (13) 

We obtain the baryon bias factor in any epoch from the above equation (12) 
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The bias factor in Eq. (14) is usually expressed by the parameter as in the 
reference [2] 



6(a) = 1 - 



(15) 



where c is a constant to be fixed by c = adec- 



We show the evolutions of (5dm (a) (solid line) and (5b (a) (dotted line) when 
ujde = —1-0 in Fig 2a. We use ^%iy[ = 0.23 and Qb = 0.04 in this figure. 
(5b (a) catches up with (5dm as expected. In Fig 2b, we also show the evolu- 
tions of baryon bias factors for the different cases. The solid lines correspond 
to the exact form of b{a)s for Ude = —0.8, —1.0 and —1.5 obtained from Eq. 
(14). All of the lines for the corresponding models are overlapped as one 
line. The dotted and dashed lines correspond to c = 0.004, 0.01, respectively 
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Figure 2: a) Evolutions of the perturbations of dark matter and baryon when 
^DM ~ 0.23 and = 0.04 for the cosmological constant models LOde = —1-0 
(from top to bottom), b) Evolutions of bias factors 6(a) from the exact 
solutions and the parameterization of it with c = 0.004 and 0.01 (from top 
to bottom). 
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C2 


Ch X 10^ 


-0.8 


-0.783555 


0.692420 


-1.99786 


-1.0 


-0.542664 


0.697915 


-1.69270 


-1.5 


-0.312706 


0.701161 


-1.22785 



Table 1: ujde is the eos of the dark energy, ci and C2 are the coefficients of (5dm 
given in Eq. (5). Chi is the coefficient of 6^. We use fi^j^ = 0.23, rig = 0.04, 
^de ~ 0.73 and a, = 0.01 to get these values. 

when we use the approximate parameterization of b{a) given in Eq. (15). 
The bias factor gets close to 1 after a ~ 0.1 because Jb gets very close to 
5dm from this epoch. This is model independent feature of 6(a) and thus the 
given approximation of 6(a) in Eq. (15) is well matched with all of models 
we choose. Thus, baryon bias factor might not be a good tool for probing 
the dark energy model except some exotic cases [6]. 

We summarize the coefficients of 5dm and 6^ used in Fig. la and Fig. 2 
in table 1. 
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Figure 3: a) Evolutions of growth factors Dg and Dws (from top to bottom) 
when = —0.4 and 0^ = 0.2. b) Evolutions of /(a) = Om(a)''' and / 
(from top to bottom) for the same values of wje and 0^. 

We also emphasize that these solutions in Eq. (5) and Eq. (6) are different 
from the well known approximate solution in the reference [7]. 

The solution given in the reference [7] is claimed to be a growing mode 
solution of Eq. (6). However, one is not able to separate the correct growing 
mode solution with decaying mode one from the general solution except 
LVde = —1/3 or —1 [9]. It is also easy to check from the initial conditions of 
the growing mode solution in Eq. (7). If one use the solution in Eq. (16), 
then one is not able to obtain the consistent coefficient value Cgw from the 
two initial conditions. We also need to interpret the solution in the correct 
way [8]. The behavior of the exact analytic solution is consistent with the 
numerical result in the reference [9]. Furthermore, the exact solutions are 
useful to extend to the other general models like modified theory of gravities 
and apply for the matter power spectrum, the mass function, the large scale 
statistics and so on. Also from the exact solution we are able to easily 
investigate the dependence of physical quantities on both model oj^e and 
the present matter energy density contrast [5]. This fact is easy to be 
ignored if one uses the well known fitting formulae [10]. 

For the comparison of the exact solution with the fitting formulae, we use 
the total matter density growth factor instead of the DM growth factor given 
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in Eq. (6). The commonly used approximation of growth index parameter 
Iws is given by [10] 
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The variation of 7^^; with respect to a is negligible. From this approximate 
parametrization 7^3, the growth factor 5ws(a) is obtained [10, 11] 



The above expression is slightly different from the original 5ws in [10] which 
is normalized to the growth factor today. It is known that the value of 
the growth factor and the growth index obtained from the above fitting 
formulae are very close to the correct ones when u^e is close to —1 [11]. 5y,s 
and fsw are good approximations for small uj^e and big O^. However, we 
need to investigate the validity of this fitting formula for general cases. As 
a counterexample, we depict the evolutions of the growth factor 5 (i.e. D) 
and the growth index /(a) = = Q.rn{ay for Ude = —0.4 model with 
fij^ = 0.2 in Fig. 3. Wc show the evolution of the exact growth factor 5g 
and that of the approximate one 5ws in the left panel of Fig. 3. The solid 
line and the dashed line correspond to 5g and 5yjs when uj^e = —0.4 and 
fij^ = 0.2. The error of the present value of the growth factor is about 17 
% in this model. Wc also show the evolution of the growth index fy^g and 
/ for the same model in the right panel of Fig. 3. The dashed line and the 
solid one correspond to f^s and /, respectively. The errors of the growth 
index for the entire epoch are smaller than those of the growth factor. It is 
about 5 % only at the present epoch. 

Thus, even though it is known as that the fitting formula is good for the 
cosmological concordance model, one should be careful for extending this 
formula to general models with other values of cosmological parameters. 
One may solve the physical quantities numerically without worrying about 
those case. However, it is useful and straight forward to use the exact 
solution if we are able to extend the exact solution to the general models 
like the time varying ujde [8]. We also are able to apply the exact analytic 
solution to probe the matter power spectrum, the evolution of mass function, 
the gravitational lensing and etc. 

We thanks S. Nesseris and I. Waga for useful comment. 
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